Hall Voltage with the Spin Hall Effect 
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The spin Hall effect does not generally result in a charge Hall voltage. We predict that in systems 
with inhomogeneous electron density in the direction perpendicular to main current flow, the spin 
Hall effect is instead accompanied by a Hall voltage. Unlike the ordinary Hall effect, we find 
that this Hall voltage is quadratic in the longitudinal electric field for a wide range of parameters 
accessible experimentally. We also predict spin accumulation in the bulk and sharp peaks of spin- 
Hall induced charge accumulation near the edges. Our results can be readily tested experimentally, 
and would allow the electrical measurement of the spin Hall effect in non-magnetic systems and 
without injection of spin-polarized electrons. 
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Currently, much attention is given to studies of 
the spin Hall effect, which allows to polarize elec- 
tron spins without magnetic fields and/or magnetic 
material^^^^^^^. In the spin 
Hall effect, electrically induced electron spin polarization 
accumulates near the edges of a channel and is zero in its 
central region. This effect is caused by deflection of car- 
riers moving along an applied electric field by extrinsic 3 
and/or intrinsic 5 mechanisms. In a non-magnetic homo- 
geneous system, spin accumulation is not accompanied 
by a charge Hall voltage because two spin Hall currents 
cancel each other 1 . The absence of Hall voltage leads to 
difficulties in probing the spin Hall effect, since measur- 
ing a charge accumulation is much easier than measur- 
ing a spin accumulation. Recently, the spin Hall effect 
has been observed both opticall y 15 ! 16 and electrically^. 
In the latter case, a charge accumulation has been cre- 
ated through injection of spin-polarized electrons into the 
sample^. 

In the present Letter, we predict that in a system 
with an inhomogeneous electron density profile in the 
direction perpendicular to the direction of main current 
flow, the spin Hall effect results in both spin and charge 
accumulations. The pattern of charge accumulation is 
determined by the interplay of two mechanisms. The 
first mechanism of charge accumulation is based on the 
dependence of spin-up and spin-down currents on local 
spin- up and spin-down densities. Spin currents, outgoing 
from regions with higher densities, are not fully compen- 
sated by incoming currents, therefore, a charge accumu- 
lation appears. This mechanism is primarily responsible 
for the non-zero Hall voltage. The second mechanism 
of charge accumulation is related to scattering of spin 
currents on sample boundaries which act like obstacles. 
Like in the case of Landauer resistivity dipoles 18 , this 
scattering leads to formation of local charge accumula- 
tion, which is also expected in traditionally-studied spin 
Hall systems. In addition, we show that in systems with 
inhomogeneous electron density the spin accumulation 
appears not only near the sample boundaries, but also in 
the bulk. Our approach does not involve any use of mag- 
netic materials and fields, therefore, the spin Hall effect 
can be measured electrically in completely non-magnetic 



system and without injection of spin-polarized electrons. 

To illustrate this effect, let us begin by considering 
a system having a step profile of electron density, as 
shown in Fig. [U There are several possible ways to fab- 
ricate such a system including density depletion by an 
electrode, inhomogeneous dopingi^, or variation of the 
sample height. What is important for us is that the per- 
pendicular (in y direction) spin currents are different in 
the regions with different electron density. Then, if we 
consider currents passing through the boundary separat- 
ing regions with different charge densities {n\ and 712), it 
is clear that the spin current from the region with higher 
electron density has a larger magnitude than the current 
in the reverse direction. The difference in currents im- 
plies charge transfer through the boundary and formation 
of a dipole layer. 

Let us now provide a quantitative analysis of this effect. 
We employ a two-component drift-diffusion model— »21, 
and in order to find a self-consistent solution, we supple- 
ment the drift-diffusion equations with the Poisson equa- 
tion. In our drift-diffusion calculation scheme, the inho- 
mogeneous charge density profile n(y) is defined via an 
assigned positive background density profile N(y) (such 
as the one in Fig. [1]), which, as discussed above, can 
be obtained in different ways. Assuming homogeneous 
charge and current densities in x direction and homoge- 
neous x-component of the electric field in both x and y 
directions, we can write a set of equations including only 
y and t dependencies: 




FIG. 1: (Color online) Spin Hall effect in a system with a 
step profile of the electron density in y direction, m > m. 
Spin currents through the boundary between m and 712 do 
not cancel each other, resulting in a Hall voltage. 
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where — e is the electron charge, ^|(n is the density of 
spin-up (spin-down) electrons, j^m) is the current den- 
sity, r 3 f is the spin relaxation time, crm) = en^^fi is the 
spin-up (spin-down) conductivity, \i is the mobility, D is 
the diffusion coefficient, e is the permittivity of the bulk, 
and 7 is the parameter describing deflection of spin-up 
(+) and spin-down (-) electrons. The current I x ,m) m 
x-direction is coupled to the homogeneous electric field 
Eo in the same direction as I x ,m) = en T(l)^o- The ^ ast 
term in Eq. ([2]) is responsible for the spin Hall effect. 

Equation (pQ) is the continuity relation that takes into 
account spin relaxation, Eq. ([2]) is the expression for the 
current in y direction which includes drift, diffusion and 
spin Hall effect components, and Eq. ([3]) is the Poisson 
equation. It is assumed that D, /i, r s f and 7 are equal 
for spin-up and spin-down electrons^ In our model, as 
it follows from Eq. ([2j), the spin Hall correction to spin- 
up (spin-down) current (the last term in Eq. [2j) is simply 
proportional to the local spin-up (spin-down) density. All 
information about microscopic mechanisms for the spin 
Hall effect is therefore lumped in the parameter 7. 

Combining equations (pQ) and ([2]) for different spin 
components we can get the following equations for elec- 
tron density n = + and spin density imbalance 



P = nt 



dn 



d_ 

dy 



finEy + D 



dn 
dy 



and 



dP 
~dt 



d_ 

dy 



liPEy + D 



dP 

dy 



■ jnfiEo 



P 

T sf' 



(4) 



(5) 



Analytical solution - Before solving Eqs. ([3])- ([5]) nu- 
merically, let us try to find analytical solutions in specific 
cases. This will help us in the discussion of the numer- 
ical results. An analytical steady-state solution of these 
equations can indeed be found for the case of exponen- 
tial density profile in a system which is infinite in the y 
direction. 

The structure of Eqs. ([3])- ([5]) allows us to select a so- 
lution in the form 



N(y) = Ae ay , 
P = Ce ay , 
E y — const, 



(6) 
(7) 
(8) 



where A, C and a are constants (A and a are assigned). 
This solution corresponds to constant spin polarization 
p = P/n. Substituting Eqs. ©-© into Eqs. (@j) and 
([5]) (note that the Poisson equation ([3]) is automatically 
satisfied) we obtain 



jiEyA + DaA + j/JjEqC = 0, 
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From these equations, eliminating E y , we find 
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The physical solution corresponds to the + sign in Eq. 
(pTJ). It can be easily verified that the solution given by 
Eqs. ©-([8]), ([TT]) corresponds to j y = 0. Substituting 
Eq. (Q3D into Eq. © we finally get 
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The second term on the RHS of Eq. ([12]) is the electric 
field needed to compensate the transverse current arising 
due to the spin Hall effect. If we now assume the sample 
has a finite (but large) width L, then, this term can be 
interpreted as due to charge accumulation near the edges, 
as in the ordinary Hall effect. Therefore, the Hall voltage 
can be written as 



V H ^L- 



1 + yj 1 + (2T sfl nE aY 
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LT s fiiaj 2 E$, 2T sf jfiE a < 1 
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From this equation we see that the Hall voltage is 
quadratic in Eo for small values of the parameter 
2T s f~fuEoa, and linear in Eq for large values of this pa- 
rameter. In fact, the quadratic dependence is quite un- 
usual, since in the ordinary Hall effect the Hall voltage 
is linear in the longitudinal current. The reason for this 
unusual dependence can be understood as follows. The 
charge current in the y direction, determined by the dif- 
ference of spin- up and spin-down currents, has a compo- 
nent (related to the last term in Eq. @) proportional to 
the spin density imbalance P times 7^0 • At small values 
of 2T s fyu'Eoa, the spin density imbalance is proportional 
to 7^0 itself. Therefore, the charge current and Hall volt- 
age are quadratic in Eo. At large values of 2r s f / yjj J Eoa^ 
the spin density imbalance saturates and the current de- 
pendence on i^o becomes linear. Another difference with 
respect to the ordinary Hall effect is that the polarity of 
the Hall voltage in the spin Hall effect is fixed by the 
geometry of the structure, and does not depend on the 
direction of the longitudinal current. 
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FIG. 2: (Color online) Distributions of the electron density 
n(y) and spin density imbalance P(y) = rij — nj, for a step (a) 
and exponential (b) background density profiles. The plots 
presented in the paper were obtained using the parameter 
values fi 8500cm 2 /(Vs), D 55cm 2 /s, e 12.4, r sf 
10ns, 7 = 10 -3 , Eq — lOOV/cm and the background density 
profiles: (a) N = 10 16 (1 + 6{y - L/2))cm~ 3 and (b) N = 
10 16 exp(2?//L)cm _3 , where #(..) is the step function, and 
L = 100/xm is the sample width. 



Let us now estimate the magnitude of 2r s f^jiEQa. 
Taking parameters related to experiments on GaAs 
(r sf = 10ns, 7 = 10- 3 (Ref. i), /i = 8500cm 2 /(Vs), 
E = lOOV/cm, a = 2/L, L = 100/im), we find 
2r s fjfiEoa = 3.4 • 10 -3 . Therefore, in experiments with 
GaAs, most likely, a quadratic Hall voltage dependence 
on the longitudinal electric field can be observed. 

Numerical solution - Equations ©-([5]) can be solved 
numerically for any reasonable form of N(y). We choose 
for their simplicity (and possibility to be realized in prac- 
tice) a step profile and an exponential profile. We solve 
these equations iteratively, starting with the electron 
density n(y) close to N(y) and P(y) close to zero and 
recalculating E y (y) at each time step^ Once the steady- 
state solution is obtained, the Hall voltage as a function 
of Eq is calculated as a change of the electrostatic poten- 
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FIG. 3: (Color online) Variations in the transverse charge den- 
sity induced by the longitudinal current. Here, 5n = n(Eo = 
lOOV/cm) — n(Eo = 0). The curve for the exponential profile 
has been shifted vertically by 3 • 10 10 cm -3 for clarity. The 
dashed lines corresponding to Sn = are there to guide the 
eye. 



tial across the sample. 

Fig. [2] shows distributions of the charge density and 
spin density imbalance in systems with a step (panel (a) 
of Fig. [2]) and exponential (panel (b) of Fig. [2]) back- 
ground densities. The values of parameters used for these 
particular simulations were selected to be close to experi- 
mental conditions reported in Ref. 15. However, we have 
tested the robustness of our predictions by solving Eqs. 
d3j)-([5]) for different values of parameters, and found that 
the predicted Hall voltage should be measurable under a 
wide range of experimental parameters. Quite generally, 
the self-consistent charge density n(y) is very close to the 
background density N(y). Small deviations of n{y) from 
N(y) can be observed in regions with strong gradients of 
N(y). In particular, we can notice that the step profile 
of electron density in Fig. [31(a) is smoothed out. Such a 
charge redistribution is related to the diffusion term in 
Eq. (jlj). The charge diffusion leads to the formation of 
a built-in electric field that equilibrates the charge diffu- 
sion. 

We also find that the induced spin density imbalance P 
in systems with inhomogeneous electron densities shows 
some new features, in addition to the well-known spin ac- 
cumulation near the edges. For instance, in Fig. Efa), P 
has an additional peak around y = 50/im. In Fig. [U[b), 
P is almost constant in the central region of the sample. 
In both cases, the physics of non-zero spin density im- 
balance is the same: the spin-current incoming from the 
right is stronger than the spin-current incoming from the 
left. We note that the integral spin density imbalance is 
always zero. 

At Eq = 0, the system is spin-unpolarized and there 
is no Hall voltage. When the longitudinal current is 
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FIG. 4: (Color online) Hall voltage as a function of the lon- 
gitudinal electric field Eo. 

switched on, the electron charge redistributes, and the 
associated Hall voltage appears. The change of electron 
density due to the spin Hall effect is presented in Fig. 
[U The first interesting observation is that there is a 
strong charge accumulation near the edges followed by 
a charge depletion region. Another observation is that 
the total electron density in the left region of the sam- 
ples (y < 50/im) has increased and, correspondingly, the 
total electron density in the right region has decreased. 
This change of the electron distribution can be seen in 
Fig. [3l Therefore, the left part of the samples is charged 
negatively and the right part is charged positively, as 
schematically shown in Fig. [U 

The mechanism of formation of sharp peaks of charge 
accumulation near the edges is similar to the mechanism 



of formation of Landauer resistivity dipoles 18 . From the 
point of view of spin currents, the sample edges act as ob- 
stacles which block the current flow, and lead to charge 
accumulation. The adjacent regions with the depleted 
electron density can be interpreted as screening clouds. 
We stress that this Landauer-type dipoles of charge ac- 
cumulation are quite general for spin Hall systems, and 
should thus be present also in traditionally studied struc- 
tures with a constant density profile. 

We finally plot in Fig. [4] the change of the electrostatic 
potential across the sample as a function of longitudinal 
electric field. The Hall voltage, for both density profiles, 
has a dependence on Eo which is very close to the ex- 
ponential dependence we have predicted analytically in 
Eq. [13] for small values of 2r s fjfiEoa. The fact that this 
exponential dependence appears also in the step profile, 
hints at a possible "general" property of the Hall voltage 
in spin Hall systems with inhomogeneous densities. We 
emphasize that a Hall voltage should also appear in spin 
Hall systems with a homogeneous electron density, but 
inhomogeneous 7. This corresponds to the case in which 
the spin-orbit coupling is dependent on spac o 24 i 25 . 

In conclusion, we have shown that a Hall voltage would 
appear in spin Hall systems with inhomogeneous electron 
density in the direction perpendicular to main current 
flow. The striking result is that this Hall voltage is gen- 
erally quadratic in the longitudinal electric field, unlike 
the ordinary Hall voltage which is linear in the same field. 
These results can be easily verified experimentally, and 
would simplify tremendously the measurement of the spin 
Hall effect by allowing an electrical measurement of the 
latter in non-magnetic systems, and without injection of 
spin-polarized electrons. 
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